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q ' Abstract. Coupled human balancing tasks are investigated based on both pseudo- 

neural controllers characterized by time-delayed feedback with random gain and 
natural human balancing tasks. It is shown numerically that, compared to single 
balancing tasks, balancing tasks coupled by mechanical structures exhibit enhanced 
stability against balancing errors in terms of both amplitude and velocity and also 
improve the tracking ability of the controllers. We then perform an experiment in 
which numerical pseudo-neural controllers are replaced with natural human balanc- 
ing tasks carried out using computer mice. The results reveal that the coupling 
ON ' structure generates asymmetric tracking abilities in subjects whose tracking abilities 
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are nearly symmetric in their single balancing tasks. 

Keywords: Neural Controller, Mechanical Coupling, Visuomotor Tracking, Sta- 
bility, Sensitivity 

1. Introduction. Competitive and cooperative dynamics can arise when multiple 
agents (autonomous entities) share common resources and environments. Extensive 
research has been conducted on such mutual interactions. Research in this area can 
be broadly classified into two categories: the field of mathematical ecology PQ, which 
finds group behavior models having low degrees of freedom, and the field of multi- 
robot systems [2J, which develops individual agents generating group behavior. In 
these conventional approaches, couplings between agents are constructed using shared 
resources or information. 

On the other hand, it is also reasonable to consider another type of coupling, con- 
structed of mechanical structures, for instance, connection rods placed between agents. 
Such a mechanical type of coupling appears to play an important role in clarifying the 
dynamics of agents, including humans and robots, which are in physical contact with 
one another. This type of knowledge will be applicable to developing robots working 
together or robots supporting human activities, based on physical contact. However, 
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it appears that little is known about the type of dynamics that is produced by mechan- 
ical coupling. Therefore, as a preliminary step toward solving this problem, we have 
developed the coupled inverted pendula (CIP) model with four DOF [3], in which non- 
linear behavior similar to interspecific competition in an ecosystem [T] can be directly 
generated by individual mechanical structures of agents. The CIP model proposed by 
the authors consists of a pair of independently PD-controlled inverted pendula, the 
tips of which are connected through a rigid rod. 

The primary approach used in the present study is to replace the PD controllers of 
the CIP model with human balancing tasks. It follows that the two subjects perform 
stick balancing tasks under the situation in which the tips of the sticks are mechanically 
connected to each other. This approach will allow us to determine how humans work 
together dynamically through physical contacts and will consequently provide helpful 
information for developing robots to support human activities. 

From an engineering point of view, stick balancing is basically achieved by the 
tracking control technique (See References (HE] for typical examples in engineering 
problems). However, human tracking control is known to exhibit further characteristics 
beyond engineering tracking control. For example, the human tracking control is 
known to have a relatively large time-delay caused by neural latency and parametric 
random fluctuations typical of neural control [BJ. In the field of statistical mechanics, 
it has been reported that random fluctuations arising in single human stick balancing 
tasks can be accurately modeled as an inverted pendulum with a time-delayed and 
randomly modulated feedback controller [BJ. Since such a controller does not simulate 
the structure of neural networks 0E], it should be referred to as a pseudo-neural 
control model of human balancing tasks. The most remarkable effect of the pseudo- 
neural controller is that, near stability boundaries, parametric noise can allow the 
controller to produce corrective movements on time scales shorter than the delay time 
of the controller under certain suitable conditions, exhibiting the scaling laws typical 
of the self-similarity dynamics of on-off intermittency [9]. This effect has also been 
experimentally confirmed in regard to physical human tasks, such as stick balancing 
[BjlTU] and visuomotor tracking on a computer screen [UJ. In these studies on the 
human balancing task, however, the case for two coupled humans has never been 
considered. Although common features of human control would remain even in the 
coupled case, there must be specific effects related to interactions between the two 
humans. This problem remains to be addressed. 

In the present paper, we consider the open problem of what kind of dynamics may 
arise when mechanical coupling is imposed on the balancing tasks by two humans 
whose dynamics can be characterized by the time-delayed and randomly modulated 
feedback controller. To achieve this, we propose a new model consisting of time- 
delayed and randomly modulated feedback controllers that are coupled by a mechanical 
structure equivalent to the proposed CIP model. For simplicity, we first derive a 
linearized reduced-order version of the CIP model, after which we replace the PD 
controllers with pseudo-neural controllers, as developed in the literature [B|IT0][11]. 
Using this model, we demonstrate the effects of coupling, namely, an improvement 
in stability and sensitivity of the corrective motions in the balancing tasks. We also 
perform an experiment in which we replace the numerical controllers with human 
visuomotor tracking tasks performed by subjects using computer mice. It is shown 
experimentally that the coupling structure between two subjects induces asymmetric 
sensitivity in the corrective motions of the subjects. 



EFFECTS OF MECHANICAL COUPLING ON THE DYNAMICS OF BALANCING TASKS 3 




Figure 1 . An intuitive example representing a coupled human balanc- 
ing task. 

2. Analytical Model. 

2.1. Single balancing tasks. Human stick balancing tasks and their scaling prop- 
erties have been accurately modeled as an inverted pendulum with a time-delayed 
feedback of random gain of the following form [B|: 

9 + ^9-asin9 + /3R(t)9(t-r) = (1) 

where r is a time delay representing the latency of neural reflexes in human balancing 
tasks and R(t) = 1 + is a random feedback gain, £(t) is standard Gaussian 

white noise, and v represents the strength of the noise. Note that a stick of length 
I and constant linear density is modeled as a = 3g/(2l), where g is the gravitational 
acceleration. 

Assuming that 9 ~ Ax, \ Ax\ <C 1, we obtain the linearized version of ([T]): 

Ax + 7 Ax - a Ax + (3R(t) Ax(t - r) = (2) 

which can be interpreted as an equation of motion of the relative displacement Ax := 
xrit) — xjnit), where xt is the displacement of the upper end of the stick and xu 
is the displacement of the lower end of the stick in the balancing task. It has been 
reported [TT] that xxit) and XM(t) are governed by the following equations: 

xt + IXt = ol Ax(t), &m + J&m — P R{t) Ax(t — t). (3) 

2.2. Coupled balancing tasks. In the present paper, we investigate the type of 
stability that can arise if two subjects that can be balanced independently are linked 
with a mechanical structure. Figure [T] shows an intuitive example representing this 
situation, in which each subject manipulates one of the sticks at the lower end along 
a mechanical slider with the goal of maintaining the stick in the upright position. 
Let qTi and qui (i = 1, 2) be the horizontal displacement of the upper and lower ends, 
respectively, of the zth stick. Then, the presence of the connecting rod can be described 
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Figure 2. Order reduction and linearization of the CIP model. 



by the distance I := q?i — 1t2 maintained by the fixed length of the rod. Note that 
this constant length yields the following equalities: 

q T1 = q T2 = ; q Tj q Tl = q T2 =: q T . (4) 

For a simpler description of this coupled task, we propose the following model: 

2q T + 2^q T = a Aq 1 (t ) + a Aq 2 (t) , 

qui + iqui = Ui(t,r), 

Ui(t,r) := Ri(t) A qi (t - t) (i = l,2) (5) 

where Ag» := qx — qui {i = 1,2) represents the horizontal displacement between 
the upper and lower ends, respectively, of the ith stick, as obtained in independent 
balancing, and Ri(t) = 1 + u£i(t) (i = 1,2) are independent random feedback gains, 
where £j (i = 1,2) are mutually independent standard Gaussian white noises. The 
second equation of (jSJ) can be reduced to the following relative form: 

Aft + 7 A®- -a(^A qi + Ag 2 ) +Ui(t,r) = (i = 1,2). (6) 

The proposed model (jSJ) provides a linearized reduced-order counterpart of the CIP 
model [3], as shown schematically in FigfSJ In other words, by assuming \8i\ 1, the 
gravitational restoring force proportional to sin 9i in the CIP model is approximated 
as a linear spring force with a negative coefficient a < 0. This setup is expected to 
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Figure 3. Largest Lyapunov exponent Ai for the single and coupled 
balancing tasks as a function of the feedback gain (3. 

provide the simplest model of coupled human balancing tasks. Note that the coupled 
model ([6]) coincides with the single model (j2J) if Agi = Aq 2 and Aqi = Aq 2 . In the 
following, we choose 7 = 50, a = 22, v = 0.6, and r = 0.1. For convenience, we refer 
to Ax and Aqi as balancing errors. 

The linearized model (jSJ) conserves the essential instability of the stick in the up- 
right position originally described in ([T]). Although this linear approximation becomes 
invalid as the balancing errors Ax and Aqi increase, these errors remain sufficiently 
small during the balancing task of maintaining the stick near the upright position. 

3. Effects of Coupling in the Balancing Model. 

3.1. On-off intermittency. The human stick balancing task exhibits on-off inter- 
mittency with respect to the balancing error [SHU], and the on-off intermittency oc- 
curs when the largest Lyapunov exponent Ai is slightly larger than zero [9J- Letting 
Q — (ft, (Zi, 92, Q2) be the solution of (jSJ) and letting Aq be the infinitesimal distance 
from the equilibrium point, we obtain the largest Lyapunov exponent Ai defined by [9J, 

Ai= limilog||Aqf(t)|| (7) 

t— >oo t 

where || A<jr|| represents a Euclid norm of the vector Aq. 

Figure [3] shows the largest Lyapunov exponent Ai of Ax and Agi as numerically 
calculated from (J2]) and respectively. The results for Aq 2 are omitted due to their 
similarity to the results for Aqi. We choose Ai = 5 x 10~ 5 to be slightly larger than 
zero in order to produce on-off intermittency of the balancing errors Ax in (T5]) and 
Aqi in (0), as shown in FigJH 

Figures 0(a) and [51(b) show double logarithmic plots of the power spectra of the 
balancing errors Ax and Ag 1; respectively. The results for Aq 2 are omitted due 
to their similarity to the results for Aqi. Both results have two linear slopes repre- 
senting two power law regimes, which is in good agreement with previous balancing 
experiments [SI EDI E], in which the two slopes were interpreted as a sign of on-off 
intermittency [9]. Based on these results, it can be concluded that the coupling term 
does not change the scaling law typical of the on-off intermittency of balancing tasks. 
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Figure 4. Balancing error Ax(t) of the single system (T5]) for Ai = 5 x 
10- 4 (13 = 20.306) and Aq 1 (t) of the coupled system © for Ai = 5 x 10~ 4 
(P = 21.032). 
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Figure 5. Scaling laws in the power spectra of (a) the single balancing 
task in (j2J) for Ai = 5 x 10~ 4 (/3 = 20.306) and (b) the coupled balancing 
task in © for Ai = 5 x 10" 4 ((3 = 21.032). 



3.2. Coupling-induced stability. The most remarkable effect of the coupling in (J5J) 
is that the balancing error is drastically suppressed by the coupling, as shown in Figs. 
4(a) and 4(b). In other words, the error Aqi for coupled balancing is as low as 1% 
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Figure 6. Logarithmic plot of the root mean squares of balancing er- 
ror generated by different samples of the white noise &(£). Solid lines 
represent the respective mean values of the logarithmic plots. 



of Ax of that of single balancing in terms of both maximal values and root mean 
square (RMS) values. Figure [6] shows logarithmic plots of 500 realizations of the RMS 
values. The mean RMS values, represented by the solid lines, are 1.5 x 10 3 for single 
balancing and 3.0 x 10 _1 for coupled balancing. The results for Aq 2 are omitted due 
to their similarity to those for Aq± . Figure [6] shows that the coupling term reduces 
the balancing error Aqi(t) to 0.02% of Ax for single balancing. This result implies 
that the mechanical coupling structure can improve the stability of the amplitude in 
human balancing tasks. 

In order to examine the velocities of the balancing errors, the ratio of the probability 
density p(Ag x ) divided by p(Ax) is shown in Fig 13 where Aqi and Ax are the 
velocity errors for the coupled balancing and the single balancing tasks, respectively. 
Figure [7] shows that p(Aq\) is approximately twice as high as p( Ax) near the origin. 
Since, in the linear approximation, the velocity errors Ax and A<ji are proportional 
to the respective slant angles of the sticks, the results suggest that the probability of 
maintaining constant angles is nearly doubled by the mechanical coupling structure. 

Based on these results, the mechanical coupling structure may reduce the balancing 
error in terms of both amplitude and velocity under suitable conditions. Although it 
is beyond the scope of the present study, one possible explanation of this coupling- 
induced stability might be that the coupling causes some stochastic averaging effect, 
whereby the two independent fluctuations £j(t) (i = 1, 2) cancel each other through the 
coupling constraint. This canceling process, however, occurs in a stochastic manner 



8 EFFECTS OF MECHANICAL COUPLING ON THE DYNAMICS OF BALANCING TASKS 




-4-3-2-101234 

Ax, 

Figure 7. Ratio of the respective probability densities of the velocities 
Aqi and Ax for coupled and single balancing. 

so that further stochastic analysis related to the noise-induced order [12] in the field 
of statistical physics will be helpful in explaining the coupling-induced stability. 

3.3. Coupling-induced sensitivity. It has been reported that near stability bound- 
aries, parametric noise can allow time-delayed feedback controllers to produce correc- 
tive movements on a time scale shorter than that of the delayed feedback [6J. In 
order to observe the effect of this improvement in sensitivity as a function of time, we 
consider the short-time cross-correlation coefficient (STCC) in the following manner. 
Letting x(t),y(t) be a pair of time series to be compared and letting At be the length 
of the time interval of short-time averaging, we define the STCC as follows: 

D , w + \ C(x -m x ,y -m y ;r,t) 

R(x, y; T)(t) = y - , (8) 

a x a y 

where 

C(x,y;r,t) : = (x(s)y(s + r)) [tit+A t], 

m x := (x(s)) [t>t+A t], Ox ■= ((x(s) - m^ 2 ) 1 ^ At] 

where (X(s))[ a: b] '■= (b — a)' 1 X(s)ds is the temporal average of X(s) over the time 
interval [a, b]. 

In the following, we draw a comparison between the velocities of the upper end and 
the lower end of each stick by using the STCC in ([8]). More specifically, we focus on 
R(xt, xm) for the single system and R(qr, Qmi) (i = 1, 2) for the coupled system 
([5]). Moreover, in order to evaluate the sensitivity of the corrective movements of the 
feedback controllers, we also define the first dominant peak points f x , f^, and Tq 2 in 
R(xt,xm) and R(qT,qMi) (i = 1,2). These peak points are assumed to evaluate the 
tracking abilities of the feedback controllers. 

Figure [S] shows the realization of the STCCs, R(qr, <?a/i) and R(qr, <?a/2) of the cou- 
pled system §5§ at t = 36 for At = 5. In this plot, there are single peaks at = 0.105 
and fg 2 = 0.035. Since the delay time of the controller is r = 0.1 in our calculations, on 
average, the first controller Ui(t, r) corrects the stick movement more slowly than the 
delay time, whereas the second controller U2{t, r) performs the corrections in a shorter 
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Figure 8. Short-time cross-correlation coefficient at t = 36. 

time than the delay time. Regarding f qi , f q2 as the tracking ability of the controller, 
the second controller can track the stick movement 0.105/0.035 = 3 times faster than 
the first controller, which indicates that symmetrically placed controllers u\ and «2 
having the same specifications can develop asymmetric tracking abilities over short 
time scales. 

The peak points f ± and f q . {i = 1,2) are plotted in FigJH] as a function of time t. 
The plots are numerically constructed from single realizations of numerical solutions 
of ([3]) and ([5]). The peak points randomly fluctuate over time, forming intermittent 
clusters of points. The sensitive region in which the peak point is smaller than the 
delay time is confirmed to be of sufficient length for physical observations. However, 
the effects of coupling are not clearly observable in Figj9j Therefore, the differences 
between the results (a), (b), and (c) in Figj9]in these plots are slight. 

In order to evaluate the effects of coupling on tracking ability, the probability densi- 
ties of the peak points f ± and f qi are shown in FigJTD] as averaged over 100 realizations 
over the time interval [0, 1200] of the numerical solutions of ([3]) and (J5J). The results for 
f(g 2 ) are omitted because it is quite similar to that for f(qi). It is clearly observed in 
Fig{10]that the density of coupled balancing, represented as small rectangles, produces 
a simple peak, which is 38% higher than that for single balancing, where the peak is 
placed at a time scale of f , which is 20% shorter than that for single balancing. 

Based on these results, it can be concluded that the mechanical coupling structure 
increases the probability of the occurrence of faster corrective movements and improves 
the tracking ability of the controllers. One possible interpretation could be that this 
noise-induced sensitivity is a version of the stochastic resonance [13J (SR), a random 
phenomenon modeled by a heavily damped particle moving randomly in a double 
well potential. The SR approach has been applied to single human balancing [Ti] . 
However, the SR approach to the coupled human balancing remains unclear and will 
be discussed later herein. 
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Figure 10. Probability densities of the first peak point f±(t) for the 
single case and (t) for the coupled case as constructed from 100 sam- 
ples in the time interval t e [0 : 1200]. 
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4. Experiment of Coupled Visuomotor Tracking on a Computer Screen. 

4.1. Experimental setup. We perform an experiment in which the numerical pseudo- 
neural controllers in (jSj) are replaced by natural human balancing tasks, as shown in 
Fig jnT a). In practice, the variable % in the single model ([3D is replaced with the 
movement of a mouse manipulated by a subject, and qMi and qui in © are replaced 
by those of two subjects. Each subject is presented with combinations of thick and 
thin lines on a computer screen, as shown in FigJTTTb). The thick lines represent the 
upper ends xt and qn as calculated from the numerical models dHJ) and (JSJ, while 
the thin lines represent the lower ends xm and qMi manipulated by the subjects. The 
screen resolution is 1,200 x 600 (pixels), where the range of the displacement [—3,3] 
in the numerical model maps to the horizontal range of pixels [1,1, 200] on the screen. 
The movements performed by the subjects xm and qMi are recorded and substituted 
into the numerical models ([3]) and (jSJ) with a sampling rate of 50 Hz, and the set of 
lines on the screen is animated at the same rate. 

The experimenter issues the following instructions to the subjects: 

• Each subject should collaborate with his/her partner in order to maintain the 
assigned pendulum, provided as an overhead view on the computer screen, in 
the upright position by manipulating the lower end of the pendulum, which is 
represented SIS db thin line. 

• The upper ends, represented as thick lines, are assumed to be connected by a 
rigid rod in such a way that a constant distance is maintained between the thick 
lines. 

• The experiment begins after a countdown performed by the experimenter, and 
ends after ten minutes or when any of the lines leaves the visible range. 

• The subject can abort the experiment at anytime. 

According to these instructions, after the countdown, the numerical simulation is 
started and the experiment begins. The experiment ends after ten minutes or when any 
of the lines xt, qn, and %, qMi leaves the visible range [—3,3]. The initial positions 
of the animated lines, i.e., the initial values of the models, are set to xt{0) = —0.5, 
9mi(0) = -0.6, q M2 (0) = 0.6, x T (0) = x M (0) = 0, q T (0) = -0.5, q m (0) = -0.6, 
?M2(0) = 0.6, and Qr(0) = <7mi(0) = (?M2(0) = 0. The model parameters are chosen to 
be a — (3 — 21, 7 = 50, and I = 1. 

The six subjects were healthy volunteers aged 21 to 24 years. The experiments were 
performed according to the principles of the Declaration of Helsinki and informed 
consent was obtained. The experimental protocol was approved by the Bioethics 
Committee of Utsunomiya University. The subjects were labeled A, B, C, D, E, and 
F, and three pairs (A,B), (C,D), and (E,F) are considered in the coupling task. 

4.2. Experimental results. First, we focus on subjects A and B and the pair (A,B). 
Figure fT2l shows the experimental results for single balancing performed independently 
by subjects A and B. The thick line, which corresponds to the velocity xt of the target 
(the upper end of the stick), slightly precedes the thin line, which corresponds to the 
subject's velocity xm, due to a combination of time delays from neural reflexes and 
computer processing. Moreover, for the case in which balancing is performed by both 
subjects, the power spectrum S(oj) of the balancing error Ax exhibits a power law 
behavior typical of neural controllers, which have two different power laws with an 
exponent of ~ —1/2 in the lower-frequency range [6] . 
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Figure 1 1 . Experimental setup of the coupled human balancing tasks. 
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Figure 12. Experimental velocities of the target xt (thick line) and of 
the subject's movement xm (thin line), and power spectra of the balanc- 
ing error Ax for single balancing performed independently by subjects 
A and B. 

On the other hand, Fig{13] shows the experimental results for the case in which the 
same subjects are coupled by a rigid rod of length I in the numerical model. Note that, 
as a result of the equalities (jl]), the velocities of both targets (thick lines) coincide. 
The thick line, which corresponds to the target velocity <jt, slightly precedes the thin 
line, which corresponds to the velocities c[mi and qui of the subject, and the balancing 
errors Aqi and Aq 2 are governed by a power law having an exponent of ~ —1/2 in 
the lower-frequency range. Based on these results, the mechanical coupling structure 
between the two balancing tasks maintains the time-delay and scaling-law properties 
typical of the independently performed single balancing tasks. 
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Figure 13. Experimental velocities of the target qr (thick line) and 
of the subject's movement qui (thin line), and power spectra of the 
balancing error Aqi for coupled balancing by the pair of subjects A 
(i = 1) and B (z = 2). 

4.3. Effects of coupling on stability and sensitivity. As discussed in Section 
13. 3[ the correlation time (the first dominant peak in STCC) allows us to evaluate 
the sensitivity of the corrective movements of the subjects, i.e., the tracking ability 
of the subjects. The experimental correlation times for single balancing performed 
independently by the two subjects (denoted as A and B) are shown in FigJTH as 
peak points f ±T = 0.14, 0.12 of STCC. The differences between the correlation times 
become slightly larger when the mechanical coupling structure is placed between the 
two balancing tasks, as shown in Fig{15l i.e., f± T = 0.12, 0.18. The correlation times 
obtained from different trials of the experiment and their averages are listed in Table 
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Figure 14. Experimental short-time cross- correlation coefficients R for 
single balancing performed independently by subjects A and B. 




Figure 15. Experimental short-time cross-correlation coefficients R for 
coupled balancing performed by pair (A,B). 



[TJ The averages in Table [T] show that the coupling structure increases the individual 
differences between the subjects. In other words, the single balancing tasks yield 
similar tracking abilities ((f) = 0.132, 0.136), whereas the tracking abilities become 
asymmetric ((f) = 0.128, 0.168) when coupled by the mechanical structure. Similarly, 
Table [2] lists the RMS of the balancing errors and their averages for estimating the 
stabilities of the balancing. The averages in Table [2] indicate that the coupling has the 
opposite effect on the RMS, as compared with the STCC, i.e., the coupling structure 
decreases the individual differences in the RMS of the balancing errors from (RMS) = 
3.62 and 5.34 for the single task to (RMS) = 3.02 and 2.50 for the coupled task. 

The results in Table [1] and [2] imply that the coupling structure increases the indi- 
vidual differences between the tracking abilities (STCC), while decreasing those of the 
balancing stabilities (RMS). This effect on the six subjects can be seen in FigJT6| which 
shows experimental plots of the correlation time r versus the RMS of the balancing 
errors over all trials of the single balancing tasks performed by subjects A, B, C, D, E, 
and F and the coupled balancing tasks performed by pairs (A,B), (C,D), and (F,D). 
The plots with filled triangles indicate the results for the single tasks, and the plots 
with rectangles indicate the results for the coupled tasks. The oval is centered at the 
average point of the corresponding plots and the radii in the horizontal and vertical 
directions represent the standard deviations of the correlation time and of the RMS of 
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Table 1. Correlation times at different trials for pair (A,B). 



Subject 




Index of trial 




Average 




1 


2 


3 


4 


5 


(r) 


Single: 














A 


0.12 


0.14 


0.14 


0.12 


0.14 


0.132 


B 


0.14 


0.12 


0.14 


0.14 


0.14 


0.136 


Coupled: 














A 


0.12 


0.12 


0.12 


0.12 


0.16 


0.128 


B 


0.18 


0.16 


0.16 


0.18 


0.16 


0.168 



Table 2. RMS of the balancing errors at different trials for pair (A,B). 



Subject 




Index of trial 




Average 




1 


2 


3 


4 


5 


(RMS) 


Single: 














A 


2.4 


2.9 


3.8 


5.5 


3.5 


3.62 


B 


4.9 


4.9 


5.8 


5.0 


6.1 


5.34 


Coupled: 














A 


2.9 


3.5 


2.4 


3.9 


2.4 


3.02 


B 


2.4 


2.8 


2.6 


2.7 


2.0 


2.50 



balancing errors, respectively. A comparison of the single tasks for the solid oval with 
the coupled tasks for the dashed oval reveals that the coupling structure decreases the 
individual differences in correlation time, while increasing the individual differences 
in the RMS. This result statistically confirms the coupling effect arising in the two 
subjects, as shown in Table [1] and El 

One possible explanation for these effects of coupling on the individual differences 
might be the separation of roles assumed by the subjects in the case that requires 
cooperation. However, similar asymmetrical properties can arise in the numerical 
controller Ui(t, r) in the second equation of (j^J), which has no ability to cooperate. 
Moreover, the numerical controllers provide feedback of the combination of their own 
past states and the shared state, i.e., Aqi = — qMi, i n the following form: 

u i = u i (Aq i ) ( = P Ri(t) A qi (t - t) in ©). (9) 

In contrast, the human subjects can also receive visual feedback of their partner's state 
via the shared computer screen in FigJTTj and, as a result, their controller model may 
appear to have the following structure: 

u i = u i (Aq 1 , Aq 2 ) (i = l,2). (10) 

Developing precise descriptions of such a human controller model will provide the 
first step toward exploring the coupling effects, such as the changes in the individual 
differences. For this purpose, it appears that approaches based on learning theory and 
evolutionary methods [T5|fT6] will be helpful. However, a detailed explanation of these 
effects is beyond the scope of the present paper and will be discussed elsewhere. 

5. Conclusion. Coupled human balancing tasks are investigated based on pseudo- 
neural controllers modeled using time-delayed feedback with random gain. It is shown 
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Figure 16. Experimental plots of the correlation time r versus the 
RMS of the balancing errors over all trials of the single balancing tasks 
performed by subjects A, B, C, D, E, and F and the coupled balancing 
tasks performed by pairs (A,B), (C,D), and (F,D). 

numerically that, compared with the case of single balancing tasks, the coupling struc- 
ture increases the stability against balancing errors in terms of both amplitude and 
velocity and improves the tracking ability of the controllers. We then perform an 
experiment in which the pseudo-neural controller in the numerical model is replaced 
with natural human balancing tasks carried out using computer mice. The results re- 
veal that the mechanical coupling structure increases the individual differences in the 
tracking abilities between the subjects, while decreasing the individual differences in 
the stabilities of the balancing errors. The proposed model and experimental method 
are expected to provide the simplest means by which to understand the cooperative 
behavior between humans sharing mechanical contacts and to provide new insight 
into welfare engineering and related fields. In the future, we intend to perform model 
identification of human visuomotor tracking tasks in order to characterize the trade- 
off mechanism between tracking abilities and stabilities produced by coupled human 
subjects. 
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